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Abstract 

Spin-charge coupling is studied for a strongly confined two-dimensional hole gas subject to a 
perpendicular magnetic field. The study is based on spin-charge coupled drift-diffusion equations 
derived from quantum- kinetic equations in an exact manner. The spin-orbit interaction induces 
an extra out-of-plane spin polarization. This contribution exhibits a persistent oscillatory pattern 
in the strong-coupling regime. 
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INTRODUCTION 



Recently, the study of spin-polarized transport in semiconductors has received much at- 
tention because of its potential applications in the field of semiconductor spintronics. Many 
authors have focused on spin-orbit interaction (SOI) that allows for purely electric manipu- 
lation of spin polarization in semiconductors. Beside this useful feature, SOI brings also into 
play the undesired spin relaxation due to the coupling between the momentum of charge 
carriers and their spin (c.f., for instance, Ref. |l|). Owing to this inhomogeneous broaden- 
ing, each elastic and inelastic scattering mechanism opens up a spin dephasing channel 2|. 
The character of spin relaxation is quite different in systems with weak and strong SOI 3|. 
In the latter case, the magnetization can oscillate even in the absence of external fields. In 
contrast, for weakly spin-orbit coupled systems, the spin polarization decays exponentially 
unless it is permanently stimulated by external fields. 

The decay of spin polarization seems to be unavoidable because of the non- conservation of 
the total spin. Nevertheless, a special persistent spin-precession pattern has been identified 
recently The infinite spin lifetime of this persistent spin helix occurs in a combined 
Rashba-Dresselhaus model at a certain wave vector that gives rise to a special spin rotation 
symmetry. Furthermore, oscillations of the nonequilibrium spin density in real space, which 
g induced by the Rashba SOI. have been reported .n a number of recent papers 
|7|. These results on robust spin oscillations certainly encourage further experimental and 
theoretical studies of long-lived spin coherence states [8] in semiconductors with SOI. 

In this paper, we focus on a strongly confined two-dimensional hole gas (2DIIG) and 
study the mutual influence of SOI and a perpendicular external magnetic field. It is well 
known that a quantizing perpendicular magnetic field appreciably changes the transport 
properties of a two-dimensional electron gas (2DEG). The quantized energy spectrum man- 
ifests in Shubnikov-de Haas oscillations of the resistivity and may lead to the quantum-Hall 
effect. Due to the SOI-induced splitting and crossing of Landau levels, a beating pattern 
arises in Shubnikov-de Haas oscillations j^, which is used to determine the SOI strength 
from the measured magnetoresistivity. Similar quantum oscillations have been identified in 
the spin- relaxation rate 10 1. Other studies Q, [l^ deal with the combined effects of 
Rashba and Dresselhaus SOI on the magnetotransport in a 2DEG. Unfortunately, compara- 
ble investigations of a 2DHG are limited although the SOI is much stronger in such systems. 
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We mention the analysis of transport equations for the 2DHG at zero magnetic field 



15| , and the treatment 



the study of spin dephasing in p-type semiconductor quantum wells 
of the spin-Hall effect [l^. 

Our work is aimed to study the spin-charge coupled motion of holes in narrow quan- 
tum wells subject to a perpendicular magnetic field. Based on a rigorous density-matrix 
approach, spin-charge coupled drift-diffusion equations are derived for the 2DHG. In order 
to focus on general physical properties of the SOI in semiconductors, we adopt the simple 
cubic Rashba model that has been used in the literature m, 120, l2l| to simulate the 

SOI in a 2DIIG. This model has the striking peculiarity that there is no coupling between 
the spin and charge components of the density matrix. One should contrast this finding 
with the linear Rashba model, which is used to study effects of SOI in a 2DEG. In this 
model, the SOI leads to a coupling between spin and charge degrees of freedom. For a 
2DHG such a coupling is exclusively induced by external fields. Here, we treat a magnetic 
field applied perpendicular to the layer. Due to this field, the charge density and out-of- 
plane spin polarization couple to each other in the 2DHG. Consequently, an inhomogeneous 
spin polarization induces charge gradients, which are accompanied by an induced internal 
electric field calculated via Poisson's equation. The most interesting feature of our approach 
is, however, the observation that the character of the magnetic-field-induced spin-charge 
coupling differs qualitatively in the weak and strong coupling regime. For weak SOI, the 
dephasing time becomes much larger than the momentum-relaxation time so that the domi- 
nating mechanism is spin diffusion. In this regime, the field-induced magnetization exhibits 
only a smooth exponential dependence on spatial coordinates. Conversely, for strong SOI, 
the ballistic spin-transport regime is established, in which oscillations of the out-of-plane 
magnetization can occur. An experimental verification of this prediction would facilitate 
the technological exploitation of these long-lived spin states for the fabrication of logical 
gates. 

BASIC THEORY 

We treat coupled spin-charge excitations on the basis of an effective-mass Hamiltonian, 



which refers to the heavy-ho 



adopted in the literature 
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thin p-type quantum wells and which has been 



2l| as an acceptable simple approximation. Our 
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model includes short-range spin-independent elastic scattering on impurities and a constant 
perpendicular magnetic field B, from which only the Zeeman splitting is considered. The 
related heavy-hole Hamiltonian of the cubic Rashba model has the second-quantized form 

fc,A fc,A,A' k,k' A 

where a}^^^ {dkx) denote the creation (annihilation) operators with in-plane quasi-momentum 
k = {kx,ky,0) and spin A. In Eq. (Q, we introduced the Fermi energy ep, the vector of 
Pauli matrices a, and the strength u of the 'white-noise' elastic impurity scattering, which 
gives rise to the momentum relaxation time r. The heavy-hole band is described by the 
dispersion relation Sk = ^^fc^/(2m). The coupling of spin states as described by 



''{ki-ki),^{ki + e),hu, , (2) 
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is due to the Zeeman splitting hujc = 9*1^33/2 and the SOI, the strength of which is 
denoted by a. In Eq. ([2]), we have k± = k^ ± iky, k^ = kcos{ip), ky = ksm{{p), and 
hhJk = ak^. Within the Born approximation with respect to elastic impurity scattering, the 
four components (/, /) = (XIa /a > XIa A' /a'^aa') of the spin-density matrix fy satisfy the 



following Laplace-transformed quantum-kinetic equations 22 



23| 



sf--{K- k)f + iQ^{k) .f = \f-f) + /o, (3) 
m T 

sf + 2(cDfc X /) - -(k • k)f + iQ^{k)f 

m 

--(f-f) + ~7n^.-^^7 + fl (4) 

r TdEk T dsk 

in which the SOI-dependent vector 

huj/^i^k) = 3oi \_{ky — k'^)Ky — 2kxkyKx, {kl. — ky)Kx " ^k^kyKy, O] (5) 

couples the spin and charge degrees of freedom to each other. The wave vector k, refers to the 
center-of-mass motion and disappears in models that refer to homogeneous spin and charge 
distributions. Initial charge and spin densities are denoted by /o = n and /o, respectively. 
The cross line over fc-dependent functions indicates an integration over the polar angle if of 
the in-plane vector k. s denotes the variable of the Laplace transformation and takes over 
the role of the time parameter t. 



By treating the kinetic Eqs. ([3]) and (jl]) in the long- wavelength limit, coupled spin- 
charge drift- diffusion equations are derived for the angle-averaged spin-density matrix (/, /). 



The method has already been applied to a 2DEG without any external fields [23|. In this 
approach, it is assumed that carriers quickly reestablish thermal equilibrium. This fact 
justifies the ansatz f{ek, k \ s) = n{ek)F{K, \ s), where n{ek) denotes the Fermi distribution 
function. Expanding the solution of Eqs. ([3]) and (jlj) up to second order in k and calculating 
the integral over the angle we obtain our main theoretical result namely the following 
spin-charge coupled drift- diffusion equations 

{s + Do^')7-T,K^J, = n, (6) 
(s + — + D,k')J^ + ToJ=U, (7) 

Tsz 

{alsT + 2Q\2sT + 1))/,. + D^TK^J^ - 2aoUJ,T{l + Dtk^)!^ = {a^ + 2n^)Tf,^o, (8) 

{a^osr + 2n\2sT + 1))/^ + D^rn'^Jy + 2aoUJ,T{l + Dtk^)J.^ = (a^ + 21]2)r/,,o, (9) 

where we used the abbreviations do = sr + 1 and Q = uj^t. The A;-dependent coefficients in 
this set of equations have the form 

D„ = 4. = 24^^n^^f^t!5^j^. (10) 

— = ' ^^ = ^7^-77^' ^0 = , 11 

4+24a|^]^+32^^ ~_ An' - 3ag 

where we introduced the diffusion coefficient D = the Bohr magneton /i^, and mag- 

netic susceptibility x- The Eqs. (El) to ([9]) completely decouple in the absence of the external 
magnetic field, when = Tq = Uc = 0. This is a peculiarity of the cubic Rashba model. 
When a perpendicular magnetic field is applied to the 2DHG, a steady-state out-of-plane 
spin polarization arises 

/f = = ^, (13) 

which couples to the charge density. For a 2DEG the situation is different. In this case. 



the out-of-plane spin polarization couples to the in-plane spin components 
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The most 



surprising feature of our solution exhibits the spin-diffusion coefficient Dz in Eq. (ITTI) . the 



form of which agrees with a recently published result 



|23| derived by an alternative ap- 



proach. This particular diffusion coefficient becomes negative for strong SOI (f2 > (7o/vl2) 



indicating an instability of the spin system. In this regime, spin diffusion has the tendency 
to strengthen initial spin fluctuations. The competition between this self-strengthening and 
spin relaxation leads to undamped spin oscillations that are characteristic for ballistic spin 
transport. Such spin oscillations result from the coupling between the charge density and the 
out-of-plane spin polarization expressed by Eqs. and (jTj). What is interesting is that this 
unusual result for Dz can only be obtained by taking into account the off-diagonal elements 
of the density matrix. (In fact, neglecting and fy in Eq. (jl]), we obtain simply Dz = D). 
Therefore, the oscillations in the strong SOI regime have a pure quantum-mechanical origin 
that is manifested in the quasi-classical Eqs. (Ej) and ([7]). Strictly speaking, this result arises 
beyond the applicability of the drift-diffusion approach {23! . 24 1. 



The time dependence of the in-p 



governed by characteristic poles 
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ane spin polarization as described by Eqs. ([HD and ^ is 



25| that are calculated from alsr + 2Vt^{2sT + 1) = 0. 



Let us treat the strong-coupling regime f2 ^ 1 for the in-plane spin polarization that is 
determined by poles at st = —3/4 ± 2iQ. Performing the inverse Laplace and Fourier 
transformations, we obtain for the spectral spin polarization the result 



fxir, k\t) = exp 



fy{r, k\t) = exp 



16Dt 
IQDt 



3t 
47 

3t 
47 



cos(2u;fct) 
t/r 



U - ^ sm{2uJkt)fyo /(327rZ}), (14) 



fyo + 



2n 



sin(2^fct)/,o}/(327rD), (15) 



t/r ' 2n 

which describes damped oscillations of an initially at r = injected spin packet. The 
external magnetic field couples initial nonvanishing in-plane spin components to each other. 
A spot like initial in-plane spin polarization could be produced in experiment by a short 
laser pule. The evolution of this initial inhomogeneous spin distribution is described by 
Eqs. (14) and (15). 



SPIN POLARIZATION FOR A STRIPE GEOMETRY 

In this Section, the magnetic-field-induced coupling between the charge distribution / 
and the out-of-plane spin polarization in a 2DIIG is treated in more detail for a stripe of 
width 2L oriented along the x axis. To this end, the steady-state solution (s = 0, uo = 1) 
of Eqs. ([6]) and (JTj) is transformed back to the representation in spatial coordinates x and 
y. Due to the considered stripe geometry, the densities are independent of x. The variation 
of the charge density f{k,y) induces a self-consistent internal electric field Ey{k,y) that is 
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calculated from the Poisson equation. This internal in-plane electric field is a by-product of 
the spin-charge coupling. Its reaction on the spin is accounted for by drift terms in Eqs. ([6]) 
and ([7]). Its phenomenological consideration in Eq. for the carrier density is ruled by 
the concept of effective chemical potential [22]. Motivated by studies of electric-field effects 
on spin transport, we introduce a similar contribution in Eq. ([7]) for the out-of-plane spin 
polarization. Putting all together, the following set of coupled equations for spin-charge 
excitations are obtained 

Dik)j\k, y) - fxEyik, y)J{k, y) - V,{k)j[{k, y) = 0, (16) 

D,ik)Jlik,y)-fiEyik,y)j[ik,y) - -±^f^{k,y) -To{k)J{k,y) = 0, (17) 

E'y{k,y) = ^(Jik,y)-nik)), (18) 

in which /i = er/m denotes the mobility and e is the dielectric constant. Primes indicate 
derivatives with respect to y. We derive an analytical solution of these equations by calculat- 
ing the lowest-order contributions in the induced electric field Ey. Within this perturbational 
schema, we make the ansatz f = n + Af and = /i^"* + A/^, where the corrections result 
from the spin-charge coupling A/, Af\ ~ Ey. In addition, hard-wall boundary conditions 
Af^{±L) = and the existence of interface charges Ey{±L) = ±Eq are assumed. We obtain 
the analytic solution 

- _ AiA2-Eoro cosh(A2L) cosh(Aiy) - cosh(AiL) cosh(A2y) -Ey _ N{y) 

ATie/e N{L) ' Eo ~ N{L)' ^ ^ 

where Ai^2 are calculated from the secular equation 

In Eq. ( |T9ll . the abbreviation 
N{y) = A2 (^D,\l - cosh(A2i:) sinh(Aiy)-Ai (^D.Xj - -^^ cosh(AiL) sinh(A2i/) (21) 

was introduced. For weak magnetic fields uJcT ^ 1, we obtain the final result 

Affp eEoxHXlLl ..^r/r .f coshiy/Lp) cosh(Aiy) dnjek) , . 

fiB{l — [aiLd) ) lcosh(L/L£)) cosh(AiL) J dsk 



with 



Ai = I/VDzTsz, A2 = V47re/in/(D£) = (23) 



A 0.04 



V 
A 

V 



-0.02 




FIG. 1: SOI-induced out-of-plane spin polarization obtained from Eq. (j22p by integrating over e^. 
(indicated by < . . . >). At zero temperature, all quantities are calculated at the Fermi momentum 
kp. Parameters used in the calculation are: kp = 0.1 nm~^, ma/h? = 2 nm, £^0 = 100 V/cm, 
and n = 10^^/cm^. The thick and thin lines refer to weak (0 = 0.173, r = 0.05 ps) and strong 
(0 = 0.34, T = 0.1 ps) spin-orbit coupling, respectively. 

where denotes the Debye screening length. Again, we meet a peculiarity of the cubic 
Rashba model for a 2DHG. The final integral over the energy Ek is easily calculated at low 
temperatures. Due to the factor dn{ek) / dsk, the field-induced spin polarization is exclusively 
determined by energies at the Fermi surface for a degenerate hole gas. Therefore, the recently 
studied inhomogeneous broadening 15| due to elastic scattering is ineffective in this regime. 
However, inelastic scattering, which we disregarded in this work, may play an essential role 
for the formation of a persistent oscillatory spin pattern at strong SOI. 

The character of the solution for the out-of-plane spin polarization mainly depends on the 
strength of the SOI. In weakly coupled systems {il < l/\/T2, > 0), the spin polarization 
exhibits an exponential dependence as shown by the thick line in Fig. 1. The self-consistent 
coupling between spin and charge degrees of freedom leads to an excess magnetization at the 
boundaries of the stripe. The picture changes dramatically, when we consider the strong- 
coupling regime {Q > l/vT2, < 0). In this case, the wave number Ai becomes imaginary 
giving rise to spin-coherent oscillations. An example for this persistent spin pattern is 
shown by the thin line in Fig. 1. Despite the included elastic scattering on impurities, 
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the spin lifetime of these oscillations is infinite in the strong-coupling regime. Moreover, 
the oscillation amplitude is considerably enhanced at the resonance AiL = {2n + l)7r/2 
with n being any integer. A similar enhancement has been predicted for the SOI-induced 
zitterbewegung [26]. This observation also remind of a Fabry-Perot interferometer in optics. 
Like the finesse of the interferometer diverges for perfect reflective mirrors, the amplitude of 
the spin oscillations becomes infinite for the above mentioned particular values of the spin- 
orbit coupling and the width of the stripe. This idealized behavior indicates that beside 
elastic scattering on impurities also other spin-relaxation mechanisms have to be taken into 
account for a more realistic description of spin excitations at strong SOI. The experimental 
observation of the interesting persistent oscillatory spin structure is certainly challenging. 
It requires a spin detection set up with a high spatial resolution (the typical wavelength of 
the oscillations is of the order of 100 nm). As the magnetic field leads to a coupling between 
the out-of-plane spin polarization and the charge density, both the induced internal electric 
field and the charge density exhibit similar oscillations in the strong coupling regime. 



SUMMARY 



We studied a 2DIIG with SOI and elastic impurity scattering under the influence of a per- 
pendicular magnetic field. Applying an exact procedure, spin-charge coupled drift-diffusion 
equations were derived from quantum-kinetic equations for the spin-density matrix. The 
magnetic field mainly causes a coupling between the out-of-plane spin polarization and the 
charge density. The character of effects that result from this coupling strongly depend on 
the strength of the SOI. For weak SOI (^7 < 1/^/12), spin diffusion gives rise to an expo- 
nential decay of an initial spin polarization. In contrast, for strong SOI, the spin transport 
exhibits ballistic character so that oscillations of the magnetization can occur. This general 
conclusion was illustrated by a treatment of the spin polarization in a stripe composed of 
a 2DHG. The magnetic field induces a background magnetization that is superimposed by 
a contribution stemming from the SOI. The excess magnetization, which results from the 
spin-charge coupling, exhibits a persistent oscillatory spin pattern for systems with strong 
spin-orbit coupling. Similar standing and propagating spin oscillations with wavelength 
down to several nanometers have been treated for thin magnetic film samples js^]. The ap- 
plication of this mechanism for spin-wave logic gates depends on whether short-wavelength 
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spin oscillations can be manipulated and detected by a suitable experimental set up. 

This work was supported by the Deutsche Forschungsgemeinschaft and the Russian Foun- 
dation of Basic Research. 
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